Searching for the Origin of CP violation in Cabibbo Suppressed D-meson Decays. 



David Atwood 

Dept. of Physics and Astronomy, Iowa State University, Ames, IA 50011 



Amarjit Soni 

Theory Group, Brookhaven National Laboratory, Upton, NY 11973 

The recent evidence of relatively large direct CP violation in D° decay at LHCb suggests that CP 
studies in the D system may become an important new avenue for understanding CP just as studies 
in the B system have proven to be. The current level of CP violation could be consistent with the 
Standard Model or, perhaps, contain evidence of new physics. A clean Standard Model prediction of 
the CP violation in these decays would, of course, be important in understanding these results but 
hadronic uncertainties makes such a prediction difficult. In this paper, we make several suggestions 
to try seek the role of new physics. We propose that the hadronic enhancement needed to attribute 
the observed CP violation in D to two pseudoscalar modes may not operate for inclusive final states 
where it is likely that we will see asymmetries at the quark level expectation provided the source 
is the Standard Model. A simple way to implement this is to search for CP asymmetries in final 
states containing K and K but where the sum of their energies is less than the energy of the parent 
D. This is meant to ensure that the event belongs to an inclusive and not an exclusive sample. We 
also propose that CP asymmetries may be enhanced in modes where the tree is color suppressed. 
In particular, the final state p°p° is of special interest because it consists of charged pions only 
and, in addition, it can have C-even P-odd triple product correlations; similarly D s — > p°K + and 
p°K* + also appear interesting. We also emphasize the use of CPT constraints leading to interesting 
correlations. We then consider how isospin symmetry can provide observables which are sensitive to 
certain classes of new physics and are small in the Standard Model. In particular, we discuss using 
isospin analysis in the decays D —¥ nn, pn and pp as well as in D s — > K*n. We also consider how 
such analysis may eventually be supplemented by information about the weak phases in D° decay. 
In order to obtain this information experimentally, we consider various methods for preparing an 
initial state which is a quantum mechanical mixture of D° and D° . This may be done through the 
use of natural D° /D oscillations; observing D° mesons which arise from Bd or B s mesons which 
themselves are oscillating or from quantum correlations in D° pairs which arise from either ip" decay 
or B-meson decay. Observing CP violation in the magnitudes of decay amplitudes should be within 
the capability of experiments in the near future, however, obtaining the weak phases through the 
methods we discuss will likely require future generations of machines due to the large statistics that 
are likely to be needed. 

PACS numbers: 11.30.Er, 12.60.Cn, 13.25.Hw, 13.40.Hq 



I. INTRODUCTION 



Recent results [1| from the LHCb provide evidence for CP violation in D-meson decays, in particular, Acp{K + K ) — 
A C p(n + Tr~) = -0.82 ± 0.21 (stat) ± 0.11(syst)% giving a 3.5 sigma signal of CP violation. CDF has reported the 
modes separately obtaining A C p(K+K-) = -0.24 ±0.22 ±0.09% and A C p{tt + k-) = ±0.22±0.24±0.11%. CDF has 
also directly measured the difference previously measured by LHCb and obtained^ Acp(K + K~) — Acp{it + iT~) = 
-0.62±0.21(stat)±0.10(syst)%. A similar result oi A c p{K + K~) - A c p{-K+ = -0.87 ± 0.41(stat) ± 0.06(syst)% 
was also recently reported by BELLE at ICHEP2012 Q. Belle H also gave A C p{K+K-) = -0.32±0.21±0.09% and 
Acp(tt + tt~) = +0.55 ± 0.36 ± 0.09%. In the LHCb result, the cancellation of experimental uncertainties between the 
two modes plays an important role in the extraction of a significant signal for the difference in the two asymmetries. 
These measurements dominate the world average for the difference given by the HFAG group @: A C p{K+K-) - 
Acp(tt + tt~) = —0.678 ± 0.147%. The Belle results are particularly significant for individual modes, because the 
leptonic environment allows better detection of these two final states, so super KEK and the Super-B Factory Q 
should be able to produce more precise results in the future especially for individual modes. 

The LHCb result for the difference in the asymmetries appears to be large compared to the Standard Modcl(SM) Q 
based on early expectations as we will discuss below. The weak phase arises in the SM from the CKM matrix. The 
relevant combination of CKM elements which gives the weak phase between the tree and penguin graphs is 



\9 W \ w 5.6 x 10~ 4 
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where we have expanded this in terms of the the Wolfenstein parametrization[l^. As discussed in Section|II]we expect 
the CP asymmetry on the quark level to be roughly of this size. For specific final states, hadronic effects will alter this 
expectation appreciably especially since charm quark is so light. Thus the central value for Aqp{K + K~) — Acp(k + tt~) 
may not be inconsistent with the SM [ll| - [T5| but this is far from a proof that the SM is fully adequate and therefore 
the role of new physics cannot yet be ruled out El, EE E3 ■ There is, in fact, the intriguing possibility that various 
models for new physics (NP) are also able to contribute significantly to direct CP violation in D-decays fill E7l - [28| . 
For this reason and others, it is important to devise observables which can distinguish between SM and NP origins 
for this CP violation. 

Broadly speaking, for singly Cabibbo suppressed (SCS) decays it is useful to divide potential models of NP into 
two categories. The first is "penguin like" where the effective Hamiltonian is strictly AI = i. Generally this kind 
of contribution will result only when the NP contributes through a gluonic penguin. The second is "tree like" where 
the effective Hamiltonian contains a AI = | component. This includes models where there are extra massive scalar 
or vector bosons which enter at tree level as well as photon, Z or W penguin topologies. In principle, electroweak 
penguins could contribute in this way but in the SM such contributions are negligible. 

Depending on whether the new physics is tree like or penguin like, different kinds of studies may help identify the 
underlying mechanism. In this paper we consider strategies which would be helpful in both cases. 

In section [II] the requirements of CPT symmetry motivates us to devise a general test for SM versus NP. If the large 
observed CP asymmetry is due to SM alone, then we would expect CP violation in a more inclusive final state to 
converge to the SM quark level expectation given in eqn. [T] Using the number of kaons as a surrogate for the number 
of s-quarks, we suggest that CP violation in inclusive KK + X final states would provide a good test of this idea. We 
also suggest that hadronic matrix element enhancements mostly occur in only exclusive two-body modes (especially 
pscudoscalars) . With that in mind a simple method is suggested to experimentally identify inclusive events. 

In section IIIII we focus on finding additional two body decay modes which are likely to also show large CP asym- 
metries in the case of penguin-like NP. The key point here is to go after color - suppressed tree modes. In selecting 
potentially useful modes, we also consider which final states are more easily detected because they appear in the final 
state as charged particles only. Using these criteria modes of particular interest include D° — > p°p°, D s — > p°K^ + 

and D° — > K^°K . In addition, the final states r\r\ 1 rjrj' and rjcf) may be interesting because the "s-quark rich" 
nature of the final states even though they contain neutrals in the final state; these could be of special interest to 
upcoming SuperB Factories. We briefly discuss radiative final states which are considered in (29j . These modes may 
show CP violation in a large class of penguin NP models, however through CPT arguments we show that Aqp is 
likely to be suppressed. 

In all of the above, CP violation is observed through Acp, a difference in decay rate between D and D to a given 
final state. For this form of CP violation a strong phase is also required. It would be very useful to also be able to 
measure the weak phase directly, independently of the strong phase. In section IIVI we propose various methods to 
measure this phase. Three methods are considered: (1) using D° oscillation just as oscillation in B mesons are used 

to measure the weak phase in B° — > ipK s ; (2) using oscillation in B-mesons where B — >■ D p° (or similar final states) 

and (3) using correlations in D°D pairs. If we assume that the phase between D° and D decay to a given final state 
is of the same magnitude as the observed value of Acp in D —> irir and D —> KK then the statistics required is ~ 10 11 
mesons for the methods using oscillation in D-mesons, oscillation in D-mesons and correlations in D-pairs originating 
from B-meson decay. In the case where the D-pairs arise in a ip" factory, then ~ 10 9 D-mesons are required. More 
realistically, weak phases an order of magnitude larger than the currently observed Acp may be observable in the 
foreseeable future. This would indicate a situation where the strong phase is small and yet there is a large NP weak 
phase. 

In section [V] we consider tests for NP based on isospin which would apply to tree-like NP models. In some cases 
isospin can be used to isolate CP violation in the AI = | channel. Since the SM predicts no CP violation in this 
channel, such a signal would indicate the presence of NP. 

Tests of this form where the magnitude of a D amplitude is compared to that of a D decay amplitude are proposed 
for D 7T7T, D — > pir , D s — > nK* and D — > pp. In the case of D — > pir there are two separate amplitudes which can be 
used, one derived only from the D° — > 7r + 7r~7r° overall reaction and one which also includes input from D + — > p + ir°; 
P°tt + . In the case of D — > pp each polarization can provide a separate test. 

For the final states irir, pp and pir, we can also combine this analysis with weak phase determination in section HVl 
then additional tests are possible where the phase of a AI = | amplitude is compared to its conjugate. Again such a 
phase would be indicative of tree-like NP. 

In Section [VI] we discuss the statistical requirements for testing the SM, particularly for the determination of weak 
phases. In Section IVlII we give our summary and conclusion. 



3 



II. CPT AND FLAVOR SYMMETRY CONSIDERATIONS 
A. Rough Estimate for Quark Level Expectations 



The effective Hamiltonian for SCS charm decay can be written: 




where \ q = V cq V* q (note that A^ + A s + A/, = by CKM unitarity). The operators are 



Q\ = (qu) V -A(cq)v-A Q\ = {q a Uf3)v-A{cpq a )v-A 
Qz = E (cu)v-A(qq)v-A Qi = ^ (caUp)v-A($pqa)v-A 

q—u,d,s q—u,d,s 

Q5 = E ( m )v-A{qq)v+A Qd = Y (~Ca u l3)v+A(qi3q a )v-A 

q—u,d,s q=u,d,s 

(3) 

The first term proportional to X s — X c i is the tree contribution and the term proportional to At, is the penguin 
contribution. If we assume that there is a strong phase difference between the tree and penguin of 4> s trong then we 
can write the CP asymmetry in the quark level process c — > ddu as [3l| : 



|^4cp(c — > ddu) 

« A 2 \ 4 T]R sin ()) s tron g ~ 6.4 x 10~ 4 sin strong (4) 

Here R is a number of order 1 which depends on the Wilson coefficients. If we neglect the mass of the s-quark and 
hadronization effects and use the one loop evolution of the Wilson coefficients given in (3^, IH[ then numerically 
R w 1.2. The resultant asymmetry at the quark level is thus expected to be about the same as given in Eqn. Q]if the 
strong phase is near maximal. As discussed below CPT implies AT(c — > ddu) = — Ar(c — > ssu) where for a given 
decay A->B, AT(A ->£?) = r(A -> B) - T(A -> B) 

This quark level result need not be the same as the CP asymmetry in any given exclusive hadronic final state. 
Due to the significant hadronic uncertainties in the formation of specific final states, in order to characterize the CP 
violation in D-meson decay it is useful to consider symmetries which may be partially respected by strong interactions. 
The most obvious such symmetry is SU (3) flavor which, unfortunately is badly broken [l7j . We will discuss the use 
of the isospin subgroup of SU(3) in section [V] For the current experimental results, a more handy subgroup of 
SU(3) to consider is U-spin since it directly relates the K + K~ and final states. If this symmetry were strictly 

observed then Br(D° -> Tr+Tr") = Br(D° -> K + K~) and A C p(D° -> tt+tt") = -A C p(D° -> K + K~). Since 
Br(L>° -> 7T+7T-) = (1.397 ± 0.026) x 10" 3 while Br(D° -> K+R-) = (3.94 ± 0.07) x 10" 3 , more than a factor 
of 2 discrepancy, it is clear that U-spin is badly broken. Due to the large experimental error in the individual CP 
asymmetries, no firm conclusion with respect to U-spin applied to CP violation can be drawn though the central 
values tend to show opposite sign. 
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B. CPT Relations 



Another symmetry which must, of course, be respected, is CPT which implies that the width of the D and D 
mesons must be the same. This means that when summed over all final states of D-meson decay, the partial rate 
differences must vanish: 



E Ar Po = ° 

X 



(5) 
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FIG. 1: The unitarity graph showing the CPT identity Eqn. [6]for the quark level SCS charm decay. Cut #1 
indicated in the figure shows the case where the decay is c — > ddu with a ssu intermediate state providing the strong 
phase. Conversely, cut #2 indicated in the figure shows the case where the decay is c — > ssu with a ddu intermediate 
state providing the strong phase. The interfering tree graphs are not shown but are implied 

In detail, this means [HI, (3|| that Ar must be exchanged between the various final states. This exchange is 
caused by rescattering between at least two final states (say X\ and X 2 ) with different strong and weak phases. If X2 
rescattering into X\ provides a strong phase for X\ it will give rise to a contribution to Ar(Xi). This partial rate 
asymmetry will be exactly canceled by the contribution to Ar(X2) proportional to the strong phase produced when 
X\ rescatters into X 2 - 

At the quark level, the SM maintains CPT in SCS decays by an exchange of Ar between c — > udd and c —t uss, 
hence 

Ar(c -> ddu) = -AT(c -> ssu). (6) 

The "double penguin" unitarity graph in Fig. [TJ shows how this compensation arises where the two cuts indicate the 
two final states. Thus, cut #1 gives a final state with ddu where one of the amplitudes has an internal loop with an 
ssu final state. The magnitude of this graph is the same as that given by cut #2 giving a ssu final state with an 
intermediate ddu but the sign is opposite due to the internal loop being on the left side in this case. 

To draw conclusions concerning specific groups of decay modes, it is useful to break down Eqn. [6] according to 
quantum numbers conserved by strong interactions, since the exchange of Ar between final states can only occur 
between states which can rescattcr into each other. Such rescattering is via strong interactions so the general statement 
in Eqn. [5] can be refined to: 

AT(P,C,I,G,S) = (7) 

where P, C, /, G, 5 are the quantum numbers, parity, charge conjugation, isospin, G-parity and strangeness respec- 
tively. 

In applying this to SCS modes, where S = in the final state, we can further classify final states according to the 
number, Nk, of kaons and anti-kaons they contain. Notice that in general Nk & {0, 1, 2, 3} because Mb < 4Mjf and 
for S=0 in particular, Nk G {0, 2} therefore 

AT(PCIG, S = 0, N K = 0) = -AT(PCIG, S = 0,N K = 2) (8) 

The left and right sides of these equations should represent ddu and s~su quark content respectively since it is expected 
that c — > udd couples dominantly to Nk = while c — > uss couples dominantly to Nk — 2 so Ar(c — > udd) ~ 
AT(N K = 0) and Ar(c -> uss) - AT(N K = 2). 

Implementing Eqn. [8] directly for each combination of quantum numbers is difficult since many of the D decay 
modes contain multiple charged pions and so the quantum numbers may be difficult to determine on an event by 
event basis. 

What may be a more practical an experimental test of CP violation is to look for CP violation in the inclusive case 
summed over P, C, /, G. In this case CPT implies: 



Ar(5 = 0, N K = 0) + Ar(7r + n) = -AT{S = 0, Nk = 2) - AT(K K) 



(9) 
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FIG. 2: The current experimental results for Acp(ir + n~) and Acp(K + K~). The vertically hatched band shows 
the T<7 region for the measurement of Acp(ir + n~) — Acp{K + K~)\ the diagonally hatched band shows the 1-er 
region for the measurement of Acp(tt + t:^); the square hatched band shows the 1-er region for the measurement of 
Acp(K + K~). The dashed line indicates the U-spin prediction that Acp(tt + tt~) + Acp(K + K~ ") = 0. The black 
wedge indicates the result where CPT is maintained between K + K~ and 7r + 7r~, i.e. AT(K + K~) + Ar(7r + 7r~) = 0. 
The dot towards the lower right is the quark level expectation from Eqn. |4] with maximal strong phase if you make 
the naive assumption that Acp(K + K~) = Acp(ssu) and Acp{~k + t:~) = Acp(ddu) 

where A means that two body pseudoscalars are not included as we explain more below. CP asymmetry in both 
AT(S = 0, Nk = 0) and Ar(S f = 0, Nk = 2) should approximate the quark level CP asymmetry in ddu and ssu 
respectively. 

If CP violation in the PP final states is due to the SM then there must be some hadronic enhancement for 
those exclusive final states such as n + ir~ , K + K~ , which we would not expect to be present in the inclusive Eqn. |H1 
Recall also that for exclusive two pseudosaclar modes, in particular, there are well known reasons to expect large QCD 
corrections, e. g. chiral enhancements. It is quite unlikely that inclusive modes will receive such largish enhancements. 
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Thus, the inclusive CP asymmetry should be smaller, ~ 6 x 1CP 4 . On the other hand, if the largeness of the CP 
asymmetry in the PP final states is due to NP, then one would expect the inclusive asymmetry to be roughly of the 
same order as the exclusive. The larger statistics of the inclusive state may provide more accurate results for this 
channel leading to an important indication of the nature of the observed CP violation. 

In practice, observing a quantity like Ar(5 = 0, Nk — 2) is subject to the problem that it is not possible to catch 
every final state. Thus it is useful to rephrase the relation AY(S = 0, Nk = 2) ~ AT(c — > ssu) as 



lim x ^i AT(S = 0, N K = 2; x ) ~ Ar(c -> ssu) (10) 

where x is some CP invariant acceptance cut on the final states and / represents the cut where all events are accepted. 
In any case AT(S = 0, Nk = 2; %) is a CP violating quantity. 

Actually, given that in the sample of inclusive (Nk — 2) final states we do not want to include the exclusive KK 
mode, a simple working definition of inclusive is all those final states in which the sum of K and K energies is less 
than the energy of the parent D. 

Eqn. [8] can be broken down further with approximate symmetries allowing us to gain some understanding of the 
pattern of CP violation in exclusive decay modes. For instance, if U-spin were a good symmetry then Eqn. [8] could 
be broken down into groups of final states related by this symmetry. In particular it would follow that Ar(7r + 7r~) = 
— AT(K + K~). Of course U-spin is broken but in [ll( figure 2 they fit the experimental data including U-spin breaking 
allowed within the SM. This fit favors a solution where \Ag p(ir + ir~) I Ac p(K + K~)\ > 1 although the statistics are 
not yet good enough to draw any firm conclusion. Analogously, in [15| (see eq 34), the above ratio of asymmetries is 
predicted to be w 1.8. This suggests a pattern where the partial rate asymmetry exchange is mostly between these two 
final states, in which case we would expect A C p(tt + tt~)/A C p{K + K~) « ~Br(K + K~)/BR(ir + ir-) « -2.82 ±0.14. 

In Figure [2] the current weighted average results for Acp(it + tt~) and Acp(K + K~) are shown as la bands on a 
Acp(7r + 7r~) versus A CP (K + K~) plot (indicated by the diagonally hatched and square hatched regions respectively) 
as well as the world average for the Acp(K + K~) — Acp{7t + it~) result (vertically hatched band). The dashed line 
indicates the U-spin result Acp{K + K~) = —Acp (ir + ir~) while the black wedge indicates the result where we assume 
that CPT is maintained within PP final states by exchange between K + K~ and 7r + 7r~, i.e. that Ar(A' + A'~) + 
Ar(7r + 7r~) = 0. For comparison the circle in the lower left corner indicates the naive expectation where we assume 
that the meson asymmetry is the same as the quark level asymmetry. In particular we suppose here that the 
quark level asymmetry is given by Eqn. [4] with maximal strong phase and that Aqp(K + K~) = Acp(ssu) and 
Acp(ir + ir~) = Acp(ddu). 

In order to facilitate distinguishing SM from NP contributions, in Section ((V]) we will discuss relations which rely 
on isospin only which, unlike the more general SU(3), should be good to the level of a few percent. 



III. CANDIDATES FOR ENHANCED CP VIOLATION FOR PENGUIN-LIKE NEW PHYSICS 



Suppose that CP violation is the result of a large amplitude A interfering with a smaller amplitude a. If we 
normalize the amplitudes in units of square root of branching ratio, then Br(D — > /) ps |A| 2 while Acp(f) oc a/ A. If 
we want to observe the CP violation with a significance of N a , the number of mesons required is N = N% / (BrA 2 CP ). 
In terms of the amplitudes then, 



N N 

N = K/(BrA% P ) cx m fap cx ^ (11) 

So that generally N depends on a but is independent of A but a smaller value of A does enhance Acp; N is not 
affected because this is at the expense of the branching ratio. Going to a mode which has smaller branching ratio 
with higher asymmetry has the advantage of reducing the effects of systematic errors and other errors which are not 
statistical in nature, all other things being equal. 

If we assume that the observed CP violation in D° — > 7r + 7r~, K + K~ is due to penguin like NP it may be that larger 
signals of CP asymmetries will be present in similar decays where the SM tree contribution is suppressed. Following 
this rationale, in this section we focus on the cases where the SM tree is color suppressed. 

Color suppression in two body final states is a pattern which is often born out in B-meson decays. For example in 
decays to charm mesons the color allowed B° —> n + D~ has a branching ratio of (2.6±0.13) x 10 -3 while the analogous 

color suppressed mode B° — > ir°D has a branching ratio an order of mag nitude smaller of (2.61 ± 0.24) x 10~ 4 . A 
similar pattern obtains for related decays. In contrast, in the case of B-meson decay to two light pseudoscalar mesons, 
color suppression fails. Thus Br(B° -> tt+tt") = (5.13 ± 0.24) x 10~ 6 while Br(B° -> 7r°7r°) = (1.62 ± 0.31) x 10~ 6 
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but with PV and PP final states, for instance Br(B° -> tt+ p~ +7r"p+) = (2.3±0.23) x 10" 5 versus Br(B° -> 7r p°) = 
(2.0 ±0.5) x 10~ 6 , likewise Sr(B° -> p+p") = (2.42 ±0.31) x 10~ 5 versus Br(B° -> p°p°) = (7.3 ±2.6) x 10~ 7 ; color 
suppression of these modes seems to hold. 

It is to be expected that color suppression is less effective in D decays because of greater non-pcrturbativc effects 
and increased meson rescattering at the charm mass scale. Indeed this appears to be the case, for example Br(D° —> 
K'-tt+) = 3.89% while Br(D° -> A"°7r°) = 2.44% showing no color supression, likewise Br(D° -> K*--k + ) = 1.73% 

while Br(D° -> K V) = 2.28% again showing no color suppression. Conversely the Kp channel does seem to show 

the effect since Br(D° -> K~ p+) = 10.8% while Br(D° -> K° p°) = 1.32%. 

In spite of the unreliable evidence that color suppression is universally operative in D-meson decays, used with care 
and caution, it may provide us with a guide in searching for other decay modes for future searches for enhanced CP 
violation. In particular, modes where the SM amplitude tends to be suppressed and a possible NP penguin amplitude 
may be enhanced may be of particular interest in NP searches. In table [I] we list all the two body Cabibbo suppressed 
decay modes of D-mesons to ground state mesons. Here we use the notation ir^*' to indicate either or p^ and 
7r(*)° to indicate either 7r°, p° or uj° . Likewise we use to indicate either cj) or rf-'' (or at least the ss component of 
the latter). 

Decays of the form D a -> tt^°K^ + , D+ -> ttW+^W, D° -> K^°K { * )0 and L>° -> ttW^W, have the tree 

color suppressed in this way. Also modes of the form D° — > K^°K^ ° have an additional suppression of the tree 
contribution due to the fact that the ddss final state quark content is not the same as produced by the tree graph. 

This is born out by the smallness of the branching ratios of D° — s- K S K S and D° — > K*°K . Therefore, most 
promising from this list is D s — > n^°K^ + . 

As an example of how color suppression can work in the realm of Cabibbo Allowed D-meson decays, consider such 
decays to two vector final states. The decay D° — » K* + p~ has no suppression and its branching ratio is 10.8 ± 0.7%. 
The related color suppressed decays D° — > K p° and D° — > K cu have branching ratios 1.58 ± 0.35% and 1.1 ± 0.5% 
respectively. 

In the table, we also enumerate the cases where the mode cascades down to a final state which contains all charged 
particles (i.e. tt^ and K^), and give their branching ratios from (36j, where known. Final states with all charged 
states will generally be easier to detect, particularly at the LHCb. Based on these criteria, D° — > p°p°, D s — > p°K + 
and D s — > p°K* + are perhaps the most favorable channels to find CP violation due to penguin like new physics. 

From this point of view, the cases of D° — > p°p° is perhaps of particular interest to search for enhanced CP violation 
due to NP. An additional feature of this VV final state is the spin degree of freedom: there are three polarization 
states: transverse parallel (-<4||), transverse perpendicular (A±) and longitudinal (Ag) . Each amplitude could have 
different CP violation. Already existing measurements [36| of the polarization fractions show that the longitudinal 
mode dominates with a fraction of 67% longitudinal. Further measurements of the angular distribution using the 
methods of [37j will allow the extraction of the phases between the polarization amplitudes. 

As discussed in [H, [39| and in [40| for the analogous B decays, a qualitatively different feature of VV final states 
is that there can be P-odd triple product observables. Such observables can lead to either CP-odd or CP-even 
correlations depending on the combination of D and D decays. If the C-even combination is formed (adding the 
triple product of D and D) then the combination is CP-odd conversely the C-odd combination is CP-even. CPT 
then requires C-odd, CP violating amplitudes to be real whereas C-odd, CP violating amplitudes need a rescattering 
strong phase. 

CP-odd observables of this form that are C-even can be formed from untagged samples of D° mesons. This is an 
advantage for e + er B-factories where the initial state is self conjugate so the D° samples obtained at such machines 
could be used directly (except for the small asymmetry in the D meson production mechanism between inclusive 
B — > D + X versus B — > D + X which will have to be determined from separate studies) . In the case of LHCb this 
would not be true since the pp initial state is not self conjugate so tagging will, in any case, be necessary. 

In some tree like NP models, even if the SM tree is colored suppressed the NP contribution is not. This will tend 
to enhance the NP contribution to CP violation. In particular, if a qq pair is produced by a color neutral object 
(e.g. a Z' or higgs like boson) then the effective Hamiltonian will have a different color structure form the SM and so 
color suppression may not apply. In addition to the two families of modes mentioned above, D s — > jr^°K^ + , and 
D° — > 7rW7r(*)°, modes of the form D + — > tt^ + <jy-*> should have enhanced CP asymmetries in this scenario. 

Treelike CP violation of this form should contribute to the AI = | channel and so the isospin analysis in Section IVl 
should reveal this though with some unknown hadronization effects and thus may reveal a contradiction with the SM. 
In this scenario, the 7T7r, pn and pp systems are particularly suited since there are enough charge distributions to allow 
isospin analysis and the final state with two neutral mesons has a color suppressed tree graph which may enhance NP 
CP violation. Of these cases, the p p state also has the advantage that it leads to an observed final state with all 
charged mesons. 

Another class of final states which may be of interest in some models of new physics are those which are rich in ss, 
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in particular if they contain <f> rj and rf. The only three such two body modes which are kinematically allowed are 
D° — > fir], D° — > rj'rj and D° — > rjrj. Of course none of these modes leads to an all charged final state but the tree 
graph is color suppressed. There is some additional suppression since the quark content only couples to the uu part 
of the rf'} wave function which makes up only 20-30% of these mesons. 

Another manifestation of penguin-like NP which could lead to CP violating signals are radiative decays. At the 
quark level such decays would proceed through c — > uj which leads to modes like D° — > p°-f, D° ory, D + p + j 
and D s — > K* + ~f. Other radiative D-meson decays which would not be expected to receive large contributions from 
short distance radiative penguins are D° K°*j {Br = (3.28 ±0.35) x 10~ 4 ), D° <j>j (Br = (2.70 ±0.34) x 10~ 5 ). 
In [29| these kind of radiative decays are discussed in the context of new physics. According to their analysis, if the 
QCD dipolc c — > u transition operators: 



Qs = 9s^UL^T a G^c R Q' 8 = 9sI ^u R a^T a G^c L (12) 

then operator evolution from the NP scale to the charm scale would lead to comparable coefficients for the electro- 
magnetic dipolc c — > u transition operators: 



& = eQu^u L a^F^c R Q' 7 = e Q u ^u R a^F^c L (13) 

Even if the coefficient of Q^f 1 were much smaller than Q§ a t the high NP scale, the coefficients may be comparable 
at the charm scale. Assuming that the observed LHCb result is due to the effects of they [2£| suggest that 

the induced coefficient of Q^p would generate Aqp in D — > pj of O(10%) provided the strong phases involved were 
maximal (See however |30j). Using vector dominance one can estimate that Br(D — > pj) rn 10~ 5 perhaps making this 
mode a good test for new physics. 

It is, however, unlikely that the strong phase will be O(l) because the kinematics of the rescattering forces there to 
be an a s correction to satisfy CPT. To see this consider the unitarity diagrams in Figure [3] for the quark level process 
c — > vrf. Diagram 1 in this figure shows the interference of a NP penguin with the SM penguin (having an internal s 
quark) at lowest order. Note that while cut #1 corresponds to the 1*7 final state, cut #2 does not correspond to an 
on shell state since the ss pair must rescatter into a single photon. There is therefore no strong phase and so there 
cannot be a CP asymmetry from this diagram. Diagram 2 shows an order a s correction to diagram where there can 
be a strong phase since now cut #2 corresponds to an on shell state. This diagram therefore can give rise to a CP 
asymmetry but that asymmetry will be suppressed by a s since there is an extra loop. Model calculations carried out 
in[30( appear to bear out this situation. 



IV. WEAK PHASE DETERMINATION BY INITIAL STATE MIXING 

In this section we will consider CP violation which arises in the phase between D° and Z)° decaying to a common 
final state /. Using the interference between the two amplitudes, the weak phase can be directly measured without 
relying on the the existence of a strong rescattering phase. Thus it is useful to find new examples of CP violation as 
well as elucidating the source of CP violation in D-meson decays. 

In all cases these kinds of measurements are difficult and most likely cannot be carried out in the near future if the 
weak phase is of the same order of magnitude of the the currently observed CP asymmetry in PP final states. It is 
possible that the small direct CP asymmetry seen in D° — > PP results from a large weak phase in combination with 
a small (i.e. < O(10%)) strong phase. If this proves to be the case then a weak phase (O(10%)) may be measured 
by the various methods considered here. A phase of this magnitude in any mode would be hard to explain in the SM 
and so would be an indication of NP. 

If the CP violation in D-mesons is presumed to be from the SM source, the weak phase measurement can also tell 
us the magnitude of the penguin contribution. To see this, recall that the SM decay amplitude receives a dominant 
contribution from the tree which has no weak phase and a penguin contribution with weak phase 7. We can therefore 
write the amplitude for any SCS decay and its conjugate as: 

A = T + Pe + ' n 
A = T + Pe-'i 

(14) 



Cut#l 

Diagram 1 



Cut #2 




Diagram 2 

FIG. 3: This unitarity graph illustrates CPT conservation for the quark level process c — >• wy due to NP. Diagram 1 
shows the lowest order interference between NP and SM where cut #1 is for the cy final state and cut #2 is for a 
ssu final state. Cut #2 cannot be on shell. Diagram 2 shows an example of an order a s correction to diagram 1 

where in contrast cut #2 can be on shell. 



where T and P are generally complex numbers because they contain strong phases. 

Thus if we know \A\ and |A| from measurement of D° and D decays to the final state as well as the phase 
9 = arg(AA*) then 



||jl|-|jl|e i9 | =2|P|sin 7 (15) 

In the context of the SM, this allows us to extract the magnitude of the penguin. \P\ using the value of 7 known 
independently from the global fit to the unitarity triangle. To the extent that the SM contribution to the exclusive 
final state can be calculated from QCD, admittedly a very challenging task, this can be compared to such a measured 
value of \P\. 

In any case, a precise determination of the weak phase in some modes may also allow us to test the SM by feeding 
it into the isospin analysis discussed in Section [Vj In particular, the weak phase in the neutral D decay can allow us 
to determine the weak phase of the AI = | channel. Such a phase would be an indication of tree-like NP. 

In this section we consider three methods to accomplish a weak phase measurement in D° decays by producing 
mixed initial states of D° and D mesons. First we consider the "conventional approach" of taking advantage of the 
natural mixing of the two flavor eigenstates, D° versus D° . Since the rate of this mixing is small compared to the D° 
lifetime, we consider using the much larger relative mixing in the B® and B® mesons if they subsequently decay to 

a D° meson. Finally we consider methods using entangled D°lf states produced at tau-charm factories where the 
meson pair arises form a if)" resonance and also using (super) B-factories where the pair arises in a two or three body 
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Decay 


Suppressed 
Tree 


Charged 
Final State 


Favored 


Total 
BR (1(T 3 ) 


D s -J- lr l'^W+ 


X 


[p° -> TT+ir-]K + 

[p° -> 7T + 7r-][A'* + -> 7T + [A" S ->• TT+TT^]] 


X 
X 


2.7 ±0.05 


D a ->■ W A' W+ 




[<j> -S- A'+A'-]A' + 
[0 -> K+K~][K* + -> 7r+[A" 3 -> 7T+7T-]] 




< 0.3 




X 


7r+[^-> A + A'~] 


X 


2.65 ±0.08 


-+ Ar«+if wu 




A"+[A S -> TT+TT"] 

A-+[¥*°^ A+O 
[A"*+ tt+[A" s ->• 7r + 7r"]][A* -> A'+tt^] 




1.98 ±0.13 
2.45:51 
5.7 ±2.3 

— 






7T+[p° ^TT + TT-] 




0.81 ±0.15 




XX 


[A" s -> 7T + 7r~][A' s -> 7T + 7T~] 
[A-*° -> A- + 7T-][A' S ^TT + TT-] 

f"^* — V A'^TT + HA' — V -7T + -71- - ! 
I_f\ ? / \ /I J |-£Vg T /( /I I 

[AT* ->■ AT+7r-][A* -> tt+A'-] 


X 
X 

Y 

X 


0.085 ±0.014 

< 0.2 

<r n ^ 
V. u.oo 

.07 ±0.05 




X 


[p° ^7r+7T-][p ^TT+TT-] 


X 


1.82 ±0.10 






7T + 7r~ 




1.400 ± .026 




X 


D u ->■ # u 


X 


1.40 ±0.12 


£)0 JCW+ifW- 




A" + A"- 

[A"*+ -> 7T+[A' S -> tt+tt-IIA'- 

A" + [A"*- -> 7r-[A" s -»■ TT+TT^]] 
[A"*+ -)■ 7r+[A" s -> 7r + 7r-]][A*" -> 7r^[A s -> tt+tt"]] 




3.96 ± .08 
2.19 ±0.1 
0.78 ±0.06 



TABLE I: The singly Cabibbo suppressed decays of D mesons to two ground state are listed. Note that the 
notation n'*^ stands for 7r + or p + ; stands for tt°, p° or w°; 4>^ stands for <f> or rf 1 ' to the extent that rf'> is 
an ss state. For each group of decays, we have indicated whether the tree contribution is color suppressed with "X" 
and if it is both color and Zweig suppressed with "XX" . The instances which lead to an all charged final state are 
listed. The favored column are decays where the tree is colored suppressed and the final state has an all charged 
final state indicated by "X" . Where the branching ratios are known from [36j we have included it in the last column; 
this is the branching ratio including the subsequent decays to the final all charged state indicated. 



decay of a B meson (e.g. B — > DD or DDK ). First, however, we quickly review the oscillation formalism for neutral 
mesons: 



A. Oscillation Formalism 

Let us consider a generic neutral flavored meson X (i.e. X = K°, D°, B d or B s ). Defining the light eigenstate 
(Xl) with mass itixl and width Txl and heavy eigenstate (Xh) with mass vtixh and width Txh, we have: 

\X L ) = p X \X) + q X \X) 

\X H ) = p X \X) - q X \X) (16) 
Thus the flavor eigenstates evolve with time tx according to 

\x{t)) phys - al \X) - ^g x \x) 

m) phys = 4 W) ^9 X \X) (17) 
where the time dependent mixing coefficients g± are given by: 

gX _ e -{im X H + ^r X H)t _|_ g— (imxL + %TxL)t 0-8) 
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Let / be a final state which both X and X can decay to. If Aj is the amplitude for X — > / and is the amplitude 
for X — > / then the time dependent rates of X and J to / are: 



dry 
1 



iW) /) 



^-r(X(t) /) 

dr x 



1 

2*" 



(C? +C?)\A*\ 2 + (CX-CX) 



^-A X f 

pX f 



q x A f 



2S v Re (^ A f A f 



{ ^Af*A* 



I ^AfAf 



2S*Im i^AT^) 



(19) 



Here Am x = m XH - m XL , AT X = T X h ~ ^xl, t x = T x t, x x = Am x /T x , Vx = AT X /(2T X ), = cos(x x r x ), 
= am(x X T X ), = cosh(yxTx), = sinh(yxTx) and z x = x x - iy x . 
In the case of D mesons, both and yo < O(10~ 2 ), so we will expand obscrvablcs to first order in xd and yjj. 
In this limit, the above time dependent rate becomes: 



— r(X(f) -> /) = e"^ 
dr x 

-^-T(X(t) -+f) = e^-- 
cLt x 



\Af\ 2 



\A^ 2 



V 



rRe 
rRe 



-iz*A?*A^ 



+ 0{x\y 2 ) 
+ 0(x 2 ,y 2 ) 



(20) 



In some of the examples below we will consider the time integrated effect of oscillation. To first order in x x , y x 
this can be accomplished by replacing the decay amplitudes with "effective" decay amplitudes: 



B 



x 



x 



,x 



(21) 



The time integrated rate for a D° meson to decay to / is given, up to first order in xx and yx , by using this effective 
amplitude "without oscillation". 
Thus, for instance, 



/■OO 

/ dT(X -> /) dr = \Bf\ 2 



0(x 2 ,y 2 ) 



dT(X^f) dr = \B f \ 2 + 0(x 2 ,y 2 ) 



(22) 



B. Weak Phases from D°/D Oscillation 

As discussed in [U |42j in D° decay to a given final state one must consider both direct CP violation and indirect 
CP violation due to D° oscillation. Conversely, assuming that the oscillation parameters are known from separate 

studies, we can use oscillation extract the phase between Aj and Aj . To do this, it is necessary to observe the time 
dependence of the decays. 

From eqn. [20lif we know xd, Ud, P D and q D we see that the constant term gives the magnitudes of the amplitudes 
\A®\ and \Af |. The slope of the decay rate gives the phase between these two amplitudes. If / is self conjugate like 
n + ir~ then such a phase difference will be CP-odd. If / is not self conjugate, such as p + ir~ then the phase will be 
a combination of CP-odd and CP-even phase differences. If both xjj and yo are non-zero, then the phase can be 
determined in this way without ambiguity. If one of these is zero, then there is a two-fold ambiguity in the phase 
determination. 
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C. Weak Phases from Bg/B q Oscillation 

Another way to accomplish the measurement of the relative phase is to look at a two body decay of a neutral 
B-meson, B q for q = d, s, to a neutral D- meson where the D-meson subsequently decays to the final state /. If we 
observe this overall reaction B — > M°[D a — > f] (where M° is a self conjugate neutral meson and B is either Bd or B s ) 
as a function of the time of the B q decay then the D state involved in the second decay will generally be a mixture 
of the flavor eigenstates. 

Of course, once the D-meson is spawned, it will oscillate as described above. In the following we will assume that 
only the B-meson decay time is observed and therefore the D-meson decay time is integrated over. 

Let us denote by Tjjb the amplitude for B —> A1°D, Tp B the amplitude for B — > M°D, T DB - the amplitude for 
B M°D and Tjjg the amplitude for B AI°D and thus we can define the effective amplitudes for B and B 
cascading down to the final state / Using the formalism in eqn. [3T]we can define the effective amplitudes: 



Df — B^Tdb 



B ? T db 



■ B f T DB 
B f T -DB 



(23) 



Thus the time dependent decay rate integrated over the D-mcson decay time as a function of the B-mcson decay 
time is given by eqn. \T§\ 



dr B 
1 



T(B(t B ) -> M[D -> /]) 



2S°Re (^DfD? 



( ^DfDf 



—T(B(t B ) -> M[D ->■ /] 



2 6 



V —D B 



2S"lm 



„B* 



(24) 



From the above equation, assuming that xb, Ub and p /q arc known then the magnitudes and relative phase of 
Df and D f can be determined. Assuming that Tjj is also known, then by inverting Eqn. [23] we determine Bf and 
B^ . As in the last section, we can then invert the relations contained in Eqn. [21] to determine the magnitudes and 



-rD 



relative phases of Af and Af . In most cases the amplitudes Tjj B m ±T DB - will dominate over Tbb and Tjyg and 



since {Bf , Bf}, differs from {Af,Af }, by O(10 2 ) the phase between Df and Df will, to a good approximation, 

be the negative of the phase between Af and A^ . 

Let us consider some particular cases of the parent B —> MD decay. In the case of B^ some candidates are 

B d -> TT°lf (Br=2.61 ± 0.24 x 10~ 4 ) and B d -> p°Tf (Br=3.2 ± 0.5 x 10" 4 ). The latter is probably easier to 
observe since p° decays to it + it~ . Indeed if the final state of the D° decay is cither -k + tt~ or p°p° then the entire 
event has an all charged final state. Other decays of this type are Bd — > r\D (Br=2.02 ± 0.35 x 10~ 4 ), B d — > rj'D 

(Br=1.25 ± 0.23 x 10" 4 ) and B d -> wlf (Br=2.59 ± 0.3 x 10" 4 ). In principle the results from these modes can be 
combined (taking into account the CP of the final state). In this case we could have an aggregate branching ratio of 

q 0* 

~ 10 . We can also consider D instead of D and that can augment the effective branching ratio. 

Another choice is to consider decays such as B d — > K S D (Br=5.2 x 10~ 5 ) and related modes but these are an order 
of magnitude smaller in branching ratio due to Cabibbo suppression. 

It is also possible to start with a B s state. The analogous decays are B s — > K S D and related processes which likely 
have roughly the same branching ratios. These would include B s — y K S D and B s —> K*°D with branching ratios at 
the 10~ 4 — 10 -3 level (note the K* would need to decay to K s tt° to collaps the B s flavor wave function) and decays 
such as B s -> 0Z?° at the 10~ 5 - 10~ 4 level. 
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D. Correlations at Charm and B Factories 

Let us now consider the case of a D°TJ pair which is initially in a single, correlated, quantum state. Let us 
arbitrarily label the mesons D\ and D 2 and consider reactions where D\ — > f; D 2 —> g where / is the state of interest 
and g is an "index" decay (which needs to be a decay state of both D° and D ), for instance / = ir + ir~ and g = K s tt° 
where the weak phase of it + it~ is to be measured and it is assumed that the weak phase in K s ir° is small and known 
(i.e. just from K s ). 

In such a scenario, the initial wave function together with the observation of D 2 — > g determines the wave function 
of the D\ state as a mixture of the flavor eigenstates. In this way we are able to observe the interference of D° and 
D° decay amplitudes to / [H, HU . 

Starting with the wave function of the meson pair: 

* = a\D t > \D 2 > +a\Di > \D 2 > (25) 
where \a\ 2 + \a\ 2 = 1, the amplitude for the combined decay (D\ — > f)(D± — >• g) is therefore 

A fg = aA f A g + aA f A g 

\A fg \ 2 = \{\A f \ 2 \A g \ 2 + \A f \ 2 \A g \ 2 ) + i(M 2 - lol^dA/l 2 !^! 2 - \A f \ 2 \A g \ 2 ) + 2Re{a<fA f A* f A* g A g ). (26) 

Thus if I Af I, \Af\, A g , A g , 9 and S are known, then the phase between Af and Af can be determined. 

This equation takes into account only the entanglement of the initial state, we can also take into account the time 
integrated neutral D oscillations by integrating |j4/ 9 | 2 to first order in x,y. As above, this is equivalent to replacing 
Af t g in Eqn. [26] with the effective amplitudes Bf^ g given by Eqn. I2T1 

The conceptually simplest example is applying this at a tau-charm factory using the method of [43|, |3| ■ In this 
case, the D° pair arises from the decay of the if) (37 70) and so in the initial state, 

1 1 

a = +—= a = =. 

An evenly mixed Di state will thus arise when \A g \ « \A g \. As an example, if we take g = K s tt° then A g = —A g 
if, as the SM predicts, there is no CP violation in this pure tree decay mode except for the well understood O(10 -3 ) 
CP violation in the mixing of the K s . (CP violation in D decay to states which contain K s has been observed in the 
related decay D + — > K s n + by BELLE [45[ and has been shown to be consistent with CP violation only due to the 
well understood mixing in the neutral kaon.) 

This method may be generalized somewhat to the case where g is a three body decay such as K s tt + tt~ . Here, the 
decay amplitude is a function of the kinematic variables. In this case we can specify the kinematics by the variables E± 
being the energies of the tt^ in the rest frame of the D°. The amplitudes A g and A g are functions of these variables: 
A g {E + ,E~) and ~A g {E + ,E~). If the decay to g is CP-invariant then the relation A g (E + , E~) = A g (E~ ,E + ). If we 
assume that A g (E + , E~) and A g (E + , E~) are known from other studies then eqn. [26] can be use to find the phase 
between Af and Af. 

Another potential way to generate correlated neutral D-mcson pairs is at a B factory through decays such as 

B + D Q D K + (Br=2.10 ± 0.26 x 10 -3 ). More generally, it should be possible to adapt this analysis to the decays 

B+ -> D*°lfK + (Br = 4.7 ± 1.0 x 10~ 3 ), B+ -> D*°D*°K+ (Br = 5.3 ± 1.6 x 10" 3 ) and B+ -> D°D*°K+ (BR 
not yet known) which may increase the statistics by a factor of ~ 5. 

Observing the Dalitz plot of B + — > D°D K + decay and fitting it to a resonance+background will give a model for 
the phase of the decay amplitude as a function of the Dalitz plot variables. 

Let us take the Dalitz plot variables to be Ed and Ed , the energies of the D° and D in the B + frame respectively 
so that the decay amplitude will have the dependency A(Ed, Ed)- Let Ef be the energy in the B + frame of state / 
and Eg be the energy of state g. The wave function of the D-meson pair is therefore given in terms of Ef and E g by 
(note / -O- g between the two equations): 



a(E f ,E g ) = A(E f ,E g )/yf\A(E f ,E g )\* + \A(E g ,E f )\* 

a(E f ,E g ) = A(E g , Ef) / yJ\A(Ef , E g )\ 2 + \A(E g ,E f )\* (27) 
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This is because there is interference between the case where D° — > f with 75° — > g and D — > f with D° — > g. Using 
Dalitz plot phases in this way is similar to a method used by BaBar to find the phase 7 in the B-meson decay to D°K 
where the D meson subsequently decays to 3ir (46j . 

V. ISOSPIN DECOMPOSITIONS 

Since isospin is a very good symmetry of strong interactions, conclusions reached based on isospin alone should 
hold quite accurately in spite of some theoretical uncertainties due to hadronic interactions. In a recent application of 
isospin to D decays |l5| it is argued that although generally isospin breaking is 0(1%) the isospin breaking contribution 
to CP violation should be second order in the isospin breaking parameter. 

The SM predicts that there is no CP violation in the AI = | channel because the contribution to this channel is 
only through the tree graph c — > ddu while the QCD penguin which has the CP violating phase is pure AI = i. In 
principle the Electroweak penguin could introduce CP violation into the AI = § channel but, as discussed below, this 
amplitude is negligibly small in D decays. 

In principle, at higher order in the SM, the electro-weak penguin (EWP) graphs could also contribute to CP 
violation in the AI = | channel. We can see, however, such contributions will be very small as follows: First, one 
expects that these will be suppressed compared to the QCD penguin by a factor of a^/a s ~ 0(1%). In the analogous 
case of B physics, the EW penguins are thought to be large in part due to enhancement oc m t which is not the case 
in charm decays. For example, in the decay B — > A' + 7r~, |Ac?p| = 9% so if we assume that all of this CP violation is 
due to EWP interfering with the tree graph, then we can crudely estimate the corresponding EWP contribution to 
the asymmetry in D decay as follows: 

A EWP (n ^+-, „ (EWP(D)) (Pengmn(B)) , 

A r r> (L> — > 7T 7T I ~ ; ; -r— ; ; ~r Anp(o — > 7T K ) 

CP v ' (Tree(D))(EWP{B)) V ' 

Penguin(B) (\V cb \\Vub\m b )(\V ub \\V us \) , (fj + K -\ 

~ Tree(B) (\V ud \\V cd \)(\Vtb\\V ts \m t ) CP( ~* * ' 
Pengmn(B) 2 m b . 

~ Tree(B) |Kb| m t ] (28) 

which suggests even a smaller contribution. There are a number of channels where we can directly test the premise 
that isospin is a good symmetry in D-meson decays to two body final states. As we discuss below, the relative phases 
in decays to pp provide a test of isospin conservation. 

Turning now to the isospin decomposition of SCS D-meson decays, we proceed in analogy to previous work in 
the case of B — > hit and related processes [47}[48|]. In our expansion we will adopt a notation similar to [48}. For 
each particular final state we will denote the isospin amplitude by A^ AI j = A^. AI j + Ap. AI j which indicates the 
amplitude for a transition through an effective Hamiltonian with isospin change AI leading to a final state of type 
/ with total isospin /. The right hand side indicates the further decomposition of the given amplitude into tree 
and penguin contributions respectively. Likewise the notation A{j = A^.^ + Ap.^ where i,j € {H — 0} indicates the 
amplitude for a decay with the indicated charge distribution. The corresponding amplitudes for D decay are indicated 
by A. 

Using this notation, we find for the tttt final state: 

A 7T7T V A 7T7T 

. j" A 7T7T 1 A 7T7T 

with the analogous relations also apply for the charge conjugate amplitudes. 
This leads to the following "isospin triangle" relationships: 

—j=A^L + Aqq — A 7 ^ = = —j=A + _ + A QQ — A_ (30) 
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FIG. 4: A sketch of the isospin triangle given in cqn. [30] using the central value for the branching ratios of the three 
modes. The two vectors in the sketch are proportional to the isospin 2 and isospin amplitudes as shown where the 

phase of the isospin 2 amplitude is arbitrarily taken to be real. 



Figure [4] shows a sketch of such a triangle where we use the central values for the branching ratios involved. In the 
sketch we also show A\ v and A% n . 

2 ,U 2 ,U 

The case of the KK final state can also be expanded in a similar way. The resulting relations are: 



aKK . 


2 i' 1 


2 ' L 




Al K = 


1 A KK -+- 


1 A KK + 

2 i- 1 + 


2 5>° 


A KK _ 


1 A KK -+- 

2 i- 1 + 


2 


2 5'° 



(31) 



In this case there are three isospin amplitudes determining three decay amplitudes so we cannot construct a triangle 
relation such as Eqn. [301 

In the case of the irir final state, we can see from Figure @] that the two isospin amplitudes have a strong phase 
between them due to rescattering. In particular this illustrates the failure of color suppression in this system; for 
color suppression to be realized, the two amplitudes must cancel, thus be colinear in the complex plane and have 
magnitudes related by y/2A%\ = AY n - 

In contrast for the K S K S case, evidentally there is considerable suppression of the branching ratio. This makes sense 
on the quark level since not only is the decay color suppressed but it is also Zweig suppressed. For this to happen, the 
isospin 1 tree amplitude and the isospin tree amplitude must cancel fairly well: \A^,*i 1 + \A 1 ^ I { x w ^A^ 1 ^ Q . Since 

the dominant tree amplitude is suppressed, this suggests the possibility that Acp in D° — > K S K S could be enhanced 
compared to the K + K~ case at the expense of the total rate. 

The pp final states have the same form as tttt except that there are three polarization states which can arise from 
the decay of a scalar, An, Aj_ and Ag. There will thus be a separate set of isospin amplitudes so the analog to the 
above decomposition is: 
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APP{i) _ /iPPW 

dPP(«) _ APP(i) , /jW(') 

^-1 — — — 75 ^ 3 o ~i 7^ " 1 n 

+ \/6 5' 2 V3 

zlPPW _ /tPPW ,|PPM fqo-, 



where i € {_L, ||,^} indexes the polarization state. 

For each polarization state we therefore also have an isospin triangle relation: 

1 appV) j_ aPpV) 4pp(0_n— 1 iwW , -tpp(0 T pp(i) f^l 
7/2 +- + oo +o — u — 77| +- 00 -o ^ ' 

For each final state, an angular analysis [37j provides the magnitude of the three polarization amplitudes and 
the cosine of the phase angle between them. Thus (up to a 2 fold ambiguity) the relative phases between these 
three amplitudes can be determined. The relation (Eqn. I33[) gives the phase between the three amplitudes with the 
same charge distribution in terms of their magnitudes (again up to a 2 fold ambiguity). Combining the two kinds 
of information, we have 18 phase differences for 9 amplitudes in D-decays (Taking into account an overall phase, the 
system is overdetcrmincd by 8 degrees of freedom) which checks the validity of Eqn. 1331 Since this relation was derived 
using isospin conservation, the validity of this symmetry is thus quantified. 
For the pir final state the decomposition is: 

"« = ^-7i- 4 i". + ^ 

Ap™ _ Vj3 . p7T 1 . p7r 1 . p7r 

0+ v/§ ^ 2 V8 ^ V2 ^ x 

A fm _ 1 APT i_ _ A P^ i i/tP" i 1 API* 

A +- \/l2 i' 2 2 I. 1 2 y/6 h° 

which in turn leads to the following pentagonal isospin relationships: 

V3 Aft - V2(A% + A p Q l) = A p *_ + A p \ + 2A P £ 

v/3 74ft = V2(T; + 7C) = A p * + + A p + *_ + 2AZ (35) 
also the AI = 3/2 contribution to the 7 = 1 final state follows from the relation: 

3 Aft = V2(A P l - A%) + 2(A p + *_ - A p _\) 

3 2ft = V5« " + 2(A P :_ - A p _\) (36) 
In the case of the decay D s — > ttK* the isospin decomposition of the amplitudes is: 



1 ATT K* j_ >/2 K * 

A+0 7! i 71 1 



\/3 5 ^3 



Thus, VHq*" = A^* + v^Sf" (37) 

In this case, the two decay amplitudes depend on two isospin amplitudes so there is no isospin triangle relation as in 
the case of 7T7t and fm. 
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A. Phases in Dalitz Plots 



In the two body decays D — >• pit and D s — > K*ir the vectors decay in turn to two pseudoscalars, p — > 7T7t and 
if* — > if7r. The final states are therefore three body Dalitz decays]!^, [5(j. The same three scalar final state 
will, in general, receive contributions from a number of different pseudo two body channels. For example in the 
case of D s — > K*tt, the two charge distributions will contribute to the same three body final state, in particular 
D s -> K* + ir° -> K°tt + tt° and D s -> it*°7r+ -> 7\ 7r°7r + . Thus the A'°7r°7r + final state receives contributions from 
both the K*°tt + and K* + ir° channels. A fit to the the distribution in the Dalitz plot variables will therefore determine 
both the magnitudes of the two body amplitudes and also the relative phase between them as well as other channels 
which contribute to this final state such as K°p + . Note that the other decay of the K* in the above will not involve 
interference between these two channels, in particular D s —> K* + tt° -4 K + tt°tt° and D s — > K*°tt + — > K + tt~tt + . 

The same situation also applies to D° — > pir which leads to the final state 7r + 7r~7r°. In this case the pseudo two 
body channels p°tt°, p + ir~ and p~ir + all contribute so in fitting the Dalitz plot one obtains the magnitude and relative 
phases of each of these channels. 



B. Standard Model Tests using Isospin 

The main test of the SM origin for CP violation in hadronic D decays which can be accomplished using isospin 
analysis is to test the SM prediction that the tree graph which is the only contribution to the AI = 3/2 Hamiltonian, 
has no phase in the Wolfenstein phase convention. Thus, assuming EWP are negligible, any CP violation in phase or 
magnitude is contained in the AI = 1/2 component which receives contributions both from the tree and the penguin. 

In each system of decays there are therefore two kinds of tests which, in principle, can be performed. 

1. The magnitude of the AI = 3/2 transition amplitude is the same for the decay and its charge conjugate. 

2. The phase of the AI = 3/2 transition amplitude is the same for the decay and its charge conjugate. 

In the tttt final state, the system is sufficiently simple to allow us to cleanly extract three isospin related CP 
asymmetries. To fully characterize CP violation in this system, a fourth quantity must be determined by a phase 
measurement of the type described in Section ITVl 

Let us denote by 5{j the partial rate difference for final state / with charge distribution ij, i.e. 5{j = \A{j\ 2 — l^fjl 2 - 
Likewise for the isospin amplitudes denote: 



sLi - i^L/i 2 -i^i//i 2 

t>[AI I;AJ J] = R e (A-aI I^AJ J ~ ^AI I^AJ ,/) (38) 

Using this notation, Eqn. [29] implies that the CP violation in the D — > tttt decays can be rewritten as: 



3 

<> + o - l d p 



6 2 Z 3 2 U 3 



Since this gives each of the observed partial rate differences in terms of three different CP violating underlying isospin 
quantities, we can invert these and obtain: 



4 

C7T7T 

3*40 



SV = 261" 



V2 



_A 7r7r 
3 +0 



zo oo — 3+° 
2 +- 2 oo 



(40) 
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As discussed in [l5| the first expression for 8%]l implies, S^l = is a test of type (1) due to the evident fact that the 

2 Z 2 Z 

decay to tt + tt° is governed only by the AI = | Hamiltonian. 

The other two combinations indicate different features of CP violation in the AI = | channel which could be 
entirely due to SM physics. As discussed in [l5j . for 5][* to be non-zero requires that there are two contributions 
to this isospin channel which have different strong phases and also different weak phases. This would generally be 
expected to be the case in the SM since both tree and penguin contribute to AI = \- It could, however, happen 
that 5VI is small due to the strong phase difference between the two contributions to AI = i being small but in this 

2 U 

case the quantity 61% i „ could be non-zero due to strong and weak phase difference between the two different isospin 
channels. 

To make this clear, consider, for example, what happens if S3Z = 8VL = but <5?J. i n ^ 0. This would imply first 

"2 "2 "2 ) 2~ 

that \A\*\ = \AV^\ and \AV^\ = \A^\ but that the phase between A\* and A\* is different than the phase between 

A\ 2 and ATg resulting from a different weak phase between the two isospin channels. 

The measurement of the phase difference between either of the neutral amplitudes and their charge conjugates, i.e. 

AV[_ versus A, or Aqq versus A 00 using the methods in Section IIVI allows the complete determination of all the 

amplitudes and therefore all the CP violation in this system. This follows from relation eqn. [30] which implies that 
the three amplitudes form a triangle in the complex plane. The phase between A 7 ^ and either of the neutral modes 
is therefore determined (up to a 2 fold ambiguity) and so the phase of AVq is known. The same is also true for the 
charge conjugate amplitudes so ultimately the weak phase between A 7 ^ and ~A™ is determined (up to a four fold 
ambiguity). This then is a test of the SM of type (2). 

Note that if the phase difference with the conjugates is measured for both ir + ir~ and 7r°7r° final states, then there 
is a consistency check for the isospin relations because the isospin triangle for D° decay fixes the phase between A+l_ 
and Aqq while the the isospin triangle for D° decay fixes the phase between A^__ and Aqq . In addition, having both 
phase measurements will resolve the four fold ambiguity with respect to the orientation of the isospin triangles. Of 
course measuring the weak phase directly with the 7r°7r° final state using the methods described above will likely be 
experimentally difficult. 

The same discussion also applies to each polarization of the final state, the pp final state. Because all of the relative 
phases of the 9 D — > pp amplitudes can be measured as discussed above (and likewise for the 9 D — > pp amplitudes) , 
if one weak phase measurement is made then all of the weak phase differences are known. In principle, there are six 
possible weak phase differences (2 modes x 3 polarizations) which can be measured in D° decay to p°p° or p + p~ so 
there are multiple checks on this kind of measurement. 

For the pp final state then we have 3 type (1) tests of the SM by comparing the magnitude of each of the p + p Q 
amplitudes with their conjugates. There are two other independent tests which can be made by comparing the phase 
differences between the amplitudes with the phase differences of their conjugates. Finally, with an absolute weak 
phase determination and the isospin relationships Eqn. [32] we can have 3 type (2) tests for the weak phase difference 
between each of the p + p° polarizations and their conjugates. 

In the case of D — > pir the phase between the three D° — > pir amplitudes can be determined from analysis for the 
Dalitz plot distributions of D —> 7r + 7r~7r°. These amplitudes are therefore a part of a more general isospin analysis 
of D — > 37r as considered in [5l[ . Using the relationship eqn. [35] we can see that A p * is determined as a linear 

combination of these three amplitudes and so is determined up to an overall weak phase. Likewise we can extract the 
charge conjugate so the SM can be tested by comparison of \A p7r \ with lAs^ol- 

Furthermore, the relation eqn. [35J gives A P f 2 as a linear combination of the related charged D-mcson decays 

and Aq 7 ^ so that the phase of these two decays relative to the neutral decays can be determined up to a two fold 
ambiguity. We can thus use eqn. [36] to find the magnitude of \A P g |. Likewise we can extract the charge conjugate of 

2 ' 

the same amplitude and so SM can be tested by comparison of \A^ | with \ As A. Thus we have two tests of type (1) 

2 ,1 2 ' 

in this system. 

We can generate the corresponding type (2) tests for both of the A J = 3/2 amplitudes in D — > pir, if we know the 
weak phase difference between at least one of the neutral modes and its conjugate using the methods of Section HVl 
Since the relative phases between all the D decays are determined by the construction above, the weak phase difference 
will then be determined. The weak phases of the other two neutral cases (all of which are found in the same Dalitz 
plot) would then provide consistency checks. 

In the case of D s — > K*ir there is just a SM check of type (1). In this case, the phase between the two amplitudes 
A 7 ^ and A{\^ can be determined from the A' + 7r~7r° Dalitz plot. Thus using eqn. [37] we obtain the magnitude of 

| -A 3^ |. Again we can test the SM through verifying \A\ K \ = \A*3 |. Unlike the above cases, there is no way to 
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determine the weak phase of this amplitude because there is no neutral decay related by isospin. 

VI. GENERAL REQUIREMENTS FOR TESTING CP VIOLATION IN SCS DECAYS OF D-MESONS 

In order to form a rough estimate of the requirements to find CP violation and test the SM through the modes above 
let us assume that the CP violation in these SCS modes is generally at the same level as seen in the SCS modes (e.g. 
7T7T and KK), on the order of 0.1 — 1%. In terms of raw statistics, a sample of 10 5 — 10 7 would be required. Since the 
branching ratios of these modes is typically 10 -3 this would mean that 10 s — 10 10 D-mesons would be required; and 
probably an order of magnitude more depending on the acceptance for various decay modes. Indeed this is roughly 
true in the LHCb results [l[ based on an integrated luminosity of 0.62/6~ 1 the yield of K + K~ was 1.44 x 10 6 and 
the yield of ir + ir~ was 0.38 x 10 6 . These results point out important challenges which must be overcome to carry out 
such studies at LHCb and more generally at other facillities. 

At the LHCb it is, of course crucial to overcome the fact that the initial state is not charge conjugate. This, of 
course, is less of a problem at e + e~ colliders such as B-factories or tau-charm factories. In any case, aside from the 
requirement of raw statistics, it is necessary to identify and tag the initial D-meson and find the various final states. 
To this end, as discussed in Section Hill final states with all charged final state particles (i.e. ir^ and K^) will be 
easier to detect. 

Determining the phase through any of the methods discussed in section lIVl may require statistics somewhat beyond 
currently planned facilities. First consider using straight D oscillation with Eqn. 1201 Obviously the first requirement 
is the ability to track the time dependence of the decay to a precision < < 1/Tjjq . The relevant terms in the time 
dependence which we need to extract is the term oc r. This term, of course, is multiplied by the relative rate of 
mixing \z\ ~ 10~ 2 . Furthermore, if the weak phase is similar to the observed level of CP violation in magnitude for 
tt + tt~ /K + K ~ then we would expect arg(AA ) ~ 1%. If this is indeed the case you would need 10 9 final states in 
order to see the decay and, since the branching ratio is 10~ 3 you would therefore need ~ 10 12 D-mesons to start with. 

Using the double oscillation method, for example eqn. [23] in the B d case where y is small and p/q = e 2l/3 , the 
relevant term would be the one proportional to S% . If there was no weak phase then this would be proportional to 
the same sin 2/3 as B — > ipK s . In effect then, we would be looking for a deviation from the SM value of this coefficient 
by 0(1%) so we would expect to need ~ 10 5 decays to preform the measurement. In the case of the irir final state, 
the combined branching ratio would be 4.2 x 10~ 7 for the channel through either D°tt° or D°p°. This gives an initial 
requirement for the number of B-mesons to be ~ 2 x 10 11 . By combining a number of modes (e.g. B° — > D w, 
B° — > D p, B° — s- D 7r etc.) it may be possible to reduce this to the 10 10 — 10 11 range. 

Using the correlation method, if we take the decay B + — > K + D°D and use the index decay g = K s ir° with the 
final state / = tt + it~ and assuming we need to observe 10 5 events, then, not including acceptance, the numbers of B 
mesons needed is 3 x 10 12 . If we broaden the method to include B + K + D*°D* (assuming a total Br=l%) and 
use as an index state g = K s ir + n~~ and a target state / = p°p°, the number is reduced to 1.5 x 10 11 . 

Using correlations at a ip" factor, the number of DD pairs required using the above assumptions with index state 
g = K s tt + tt~ and a target state / = p°p° is 5 x 10 9 . 

It seems then that each of these methods requires an input of ~ 10 11 mesons if the phase is of the same order of 
magnitude as Acp for 7T7t and KK. This is probably beyond the capability of machines in the foreseeable future. 
If, however, the CP violating phase is an order of magnitude larger than Aqp (i.e. because the strong phase was 
O(10%)) then these requirements would be reduced by 2 orders of magnitude and perhaps such experiments could 
become possible at super B factories or the LHCb upgrades. 

Perhaps the cleanest environment to measure such phases would be at high luminosity charm factories where ~ 10 10 
meson pairs would be needed if the phase is 0(1%). Again if the phase were 10% this would be reduced by two orders 
of magnitude. 

VII. SUMMARY AND CONCLUSION 

D mixing is unique as its the only charge 2/3 bound system providing us with a great opportunity to search for 
new physics. In many interesting BSM scenarios enhanced mixing and also enhanced CP asymmetries arc expected; 
warped extra dimension models are a well known example. The recent discovery of direct CP violation in _D°-decays 
by the LHCb collaboration gives a huge impetus to these searches. The observed CP asymmetry of O(0.5%) is 
somewhat bigger than some estimates though it seems SM explanation is quite plausible. Hadronic uncertainties 
make precise predictions exceedingly difficult, therefore, for now, possible role of new physics cannot be ruled out. 
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More experimental information may well be pivotal in this instance This is the basic rationale behind this work leading 
us to make several suggestions. 

We suggest that the observed enhanced effects due to non-pcrturbativc physics may be most pronounced for the 
exclusive two pseudoscalar modes only, e.g. tttt and KK. For multipaticle (inclusive) final states, the quark level CP 
asymmetry of about 6 x 10 -4 may be releava nt. A simple way to implement this experimentally may be to look for 
(say) , decays of D to final states with a K and a K where the sum of their energies is less than the energy of the 
parent D. If these inclusive final states also show enhanced CP asymmetries (say at the level seen in exclusive K + K~ , 
7r + 7f~), then that would mean that it has a new physics origin otherwise it will give support to a SM explanation. 

Since the tree contribution is likely suppressed in color-suppressed final states, it is likely that CP asymmetries 
will be enhanced therein. To fecilitate experimental detection final states leading to charged its may be best to focus 
on. These twin considerations lead us to suggest D a —> p°p°, D s — > pK + and D s — > p°K* + especially interesting. 
The vector vector final states have the additional bonus that angular correlations can also be used for additional 
CP-violating observablcs. 

The importance of CPT constraints on CP violating obscrvables are emphasized and illustrated with regard to 
exlusive and inclusive and radiative modes. 

While SU(3) and Uspin symmetries seem quite badly broken in D decays, isospin likely holds quite well motivating 
us to to investigate its use especially in decays such as D — > kit, pn, pp as well as for D s — > K*tt. 

We also studied how such analysis may be augmented by information about the weak phase in D° decays. To do 
this it is necessary to study a sample of D-mesons which are in a mixed state of D° and D . Such a state may result 
from D a D oscillation or from the decay of a B or B s meson which itself is in a mixed state due to its oscillation. 
Alternatively, if a D°D pair is in an entangled state, the observation of the decay of one neutral D-meson implies 
the other D-meson is in a mixed state. Such entangled pairs may be produced in charm factories through the decay 
of ip or as the result of B-meson decays. 
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